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ABSTRACT 

Using a sample of gamma-ray burst (GRB) afterglows detected by both the X-Ray 
and the UV/Optical Telescopes (XRT and UVOT) on Swift, we modelled the spectral 
energy distributions (SEDs) to determine gas column densities and dust extinction in 
the GRB local environment. In six out of seven cases we find an X-ray absorber asso- 
ciated with the GRB host galaxy with column density (assuming solar abundances) 
ranging from (0.8 — 7.7) x 10^^ cm~^. We determine the rest-frame visual extinction 
Ay using the SMC, LMC and Galactic extinction curves to model the dust in the 
GRB host galaxy, and this ranges from Ay = 0.12 ± 0.04 to Ay = 0.65lo;o7. The 
afterglow SEDs were typically best fit by a model with an SMC extinction curve. In 
only one case was the GRB afterglow better modelled by a Galactic extinction curve, 
which has a prominent absorption feature at 2175 A. We investigate the selection ef- 
fects present in our sample and how these might distort the true distribution of Ay in 
GRB host galaxies. We estimate that GRBs with no afterglow detected blueward of 
5500 A have average rest-frame visual extinctions almost eight times those observed 
in the optically bright population of GRBs. This may help account for the ~ 1/3 of 
GRBs observed by Swift that have no afterglow detected by UVOT. 

Key words: gamma-rays: bursts - gamma-ray: observations - galaxies: ISM - dust, 
extinction 



1 INTRODUCTION 

There is now a wealth of observational evidence link- 
ing long duration gamma-ray bursts (GRBs) (prompt 
emission lasting > 2 s) with the collapse of a mas- 
sive star (coUapsar model; e.g. IWooslevI Il993l ). This in- 
cludes the underlying supernova features in the after- 
glow of some GRBs [e.g. GRB 9804 25 l|Kulkarni et al.l 
1 19981). GRB 30329 (Hiort h et all 120031 ') and GRB 060218 
I Campana et a l. 2006)], and the association b etween GRB 
host galaxies and high-mass star formation (|Tanvir et al.l 
120041 ). The 7-ray energy emission of GRBs is unaffected by 
dust and gas in the intervening interstellar medium, which 
combined with the vast energy released during the initial ex- 
plosi on, allows them to be d etected out to very high redshifts 
fe.g. lTagliaferri et al.llioosi '). These two facts give GRBs the 
potential to be highly powerful tools with which to trace 
the star formation history (SFH) in an unbiased way. Fur- 
thermore, the longer lived, lower energy afterglows light up 



their host galaxies, albeit for only a brief time (on the or- 
der of weeks) , providing invaluable insight into the chemical 
makeup of these galaxies that would otherwise be unattain- 
able in a majority of cases; certainly at high redshifts. 

To fully maximise the potential that GRBs offer as cos- 
mic probes the selection effects present in GRB studies need 
to be known, such as the stellar populations that they trace 
and the properties of their host galaxies. Infra-red observa- 
tions of long GRB host galaxies taken with Spitzer indicate 
that they are not strong starburs t galaxies, nor are they 
particularly dusty ijbe Floc'h et al.ll2CI06l ). This is in conflict 
with expectations from the coUapsar model, which predicts 
GRBs to occur in regions of active star formation that are 
heavily enshrouded by dust. This mismatch between obser- 
vations and theory could very possibly be the result of selec- 
tion effects, whereby the near infra-red (NIR) to ultra violet 
(UV) afterglow of GRBs in dusty galaxies are extinguished, 
thus reducing the chance of identifying the host. Further 
detailed analysis of GRB environments is important to de- 
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termine the range in host galaxy properties, and, fundamen- 
tally, to provide a better understanding of our observations 
and the limitations that they present. 

Due to the broadband power law spectral behaviour of 
GRB afterglows, the effects of absorbing dust and gas in the 
local environment on the GRB spectral energy distribution 
(SED) can be well identified. In the ana lysis of X-ray and 
optica l afterglow spectra for eight GRBs, Gal ama fc Wiiera 
l|200lll found evidence for high column densities of gas in the 
GRB local environment that were comparable with those ob- 
served in giant molecular clouds. Furthermore, they found 
the optical extinction to be 10-10 times smaller than ex- 
pected given the column densities. IStratta et al.l l|2004 ) ex- 
tended this sample and found the ratio of host galaxy col- 
umn density to visual extinction to be an order of magnitude 
larger that that observed in the Milky Way (MW), and also 
greater than that in the Small Magellanic Clouds (SMC) 
and Large Magellanic Clouds (LMC). Prompted by theoret- 
ical studies that indicate that the intense radiation emitted 
by a GRB should destroy small dust grains out to radii of 
around 10 pc [e.g. Fruchter, Krolik & Rhoads (2001); Perna 
& Lazzati (2002); Perna, Lazzati & Fiore (2003)], the large 
gas-to-dust ratio was taken to be evidence of the destruc- 
tion of dust in the surrounding vicinity of the burst, thus 
red ucing the visua l extin ction observed. 

IStratta et al.l (|2004l ') also found little evidence of the 
strong 2175 A Galactic absorption feature in their optical- 
NIR spectral analysis. Instead they found the SED to be 
best fit by a model in which the host galaxy has an SMC 
or starburst galaxy dust extinction law, which has no such 
absorption feature. More recently Kann, Klose & Zeh (2006) 
analysed the afterglow spectral energy distributions in the 
optical and NIR bands on a larger sample of 30 pre- Swtft 
GRBs, and also found the SMC extinction curve to provide 
a better fit to the spectra than the MW or LMC extinction 
curves. 

The simultaneous observations taken with the X-Ray 
Telescope [XRT; Burrows et al. (2005a)] and UV/Optical 
Telescope [UVO T; Roming et al. (2005)] on-board the 
Swift spacecraft l|Gehrels et al.l |2004| ) provide GRB after- 
glow SEDs without the need to extrapolate data to the same 
epoch. This capability is unique to Swift and the accurate 
broadband spectral modelling that is possible with this pro- 
vides an important data set to compare to previous multi- 
wavelength GRB samples. In this paper we model the SEDs 
of a sample of Swift GRBs with afterglows detected by both 
the XRT and UVOT to investigate the rest-frame visual ex- 
tinction and soft X-ray absorption at the GRB host galaxy. 
This provides an indication of the dust and gas content in 
the local environment of the GRB, which we compare to the 
environment of the MW, the LMC and the SMC. 

The effects of the selection bias that is introduced by 
analysing only those GRBs with X-ray and optical after- 
glows is addressed by using our results to examine the dust 
in the local environments of GRBs that have no UVOT de- 
tected afterglow. Our results have implications for the role 
that dust plays in accounting for the lack of an UV/optical 
afterglow detection in ~ 1/3 of GRBs observed by the 
UVOT. 

In section [5] we present the GRB sample and describe 
the X-ray and UV/optical data reduction and analysis and 
in section |3] we describe the models that we used to fit the 



data. We present the results from our spectral modelling 
and discuss their implications in sections |4] and [5] and also 
investigate how GRBs that have no UV/optical afterglow 
detected by UVOT may be affected by dust. Our conclusions 
are summarised in section[6l Throughout the paper temporal 
and spectral indices, a and /3 respectively, are denoted such 
that Ft oc and F^, oc and all errors are la unless 

specified otherwise. 



2 DATA REDUCTION AND ANALYSIS 

SEDs at a single-epoch were produced for a total of seven 
bursts. This epoch was chosen to be 1 hour after the onset 
of the initial prompt emission (time T) for four of the bursts 
(GRB 050318, GRB 050525, GRB 050802 and GRB 060512). 
For the remaining three GRBs (GRB 050824, GRB 051111 
and GRB 060418) the SEDs were produced at a time T-|-2hrs 
due to a lack of data at T+lhr, resulting from spacecraft 
observing constraints. 

The GRBs in our sample were chosen on the basis that 
they had afterglows detected by both the XRT and UVOT 
instruments and a spectroscopic redshift with z < 1.75. 
Bursts for which the photometry was considered to be too 
poor to provide useful constraints on the spectral fitting 
were not used. This included those with UVOT detections 
in fewer that three filters, or bursts that did not have well 
enough sampled light curves with which to obtain reliable 
multi-band photometry at a single epoch. To avoid any con- 
tamination from absorption caused by the Lya forest only 
UVOT data with a rest-frame wavelength A > 1215 A were 
used in the spectral analysis. The redshift limit was, there- 
fore, a necessary requirement to ensure that there were suf- 
ficient optical and UV data points in the afterglow SED to 
constrain the spectral fitting. At the redshift of our sample 
the afterglow spectra cover the wavelength range where the 
redshifted 2175 A extinction bump is expected to lie, and 
those GRBs closer to the redshift upper limit also probe the 
rest-frame far-ultraviolet (FUV) spectra, where the diver- 
gence between extinction curves is greatest. 



2.1 UVOT Data 

The UVOT contains three optical and three ultra-violet 
lenticular filters, which cover the wavelength range between 
1600 A and 6000 A. Photometric measurements were ex- 
tracted from the UVOT imaging data using a circular source 
extraction region with a 6" radius for the V, B and U optical 
filters and a 12" radius for the three UVOT ultra-violet fil- 
ters to remain compatible with the current effective area cal- 
ibrations^. Where possible the background rate was taken 
from an annulus with 12" inner radius and 20" outer ra- 
dius centred on the source. In the cases where there were 
nearby sources that contaminated this extraction region, the 
background was taken from a source-free region close to the 
target with radii ranging between 10 and 20". 

For each of the UVOT lenticular filters the tool 
UVOt2pha (version 1.1) was used to produce spectral files 



^ http: / /heasarc. gsfc.nasa.gov/docs/heasarc/caldb/swift/docs/uvot/ 
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Table 1. UVOT Data 



GRB 


z 


Galactic E(B-V) ^ 


Oopt 


050318 


1.44 2 


0.017 


0.94 9 


050525 


0.606 ^ 


0.095 


1.56,0.62 1" 


050802 


1.71 


0.021 


0.82 " 


050824 


0.83 5 


0.034 


0.55 12 


051111 


1.549 ^ 


0.162 


0.96 13 


060418 


1.49 


0.224 


1.25 


060512 


0.443 8 


0.014 


0.84 



Schlegel et aP lll99^) 

Berger fc MulchaevI ||2005| ') : ^iFolev et al.l ||2005|'): 



(I2005al): ^ iFvnbo et alj ll2005bl'): " iProchaska 
iDupree et al.l tOOdl ): ^ iBloom et aU (12009) 
9 IStill et al.l I2OO5I) 



Fvnbo et al.l 
l|2005l '): 



Light curve best fit by double power law l lBlustin et al.ll2006l) 
Oatcs ct al. ( 200^) 
Halpern fc MirabSl ||2005|') 
Butler et al .1 l|2006h 



Table 2. X-Ray Data 



GRB 


Data 
Mode 


Gal. Nh 
(10^0 cm-2) 


Start 

Time (s) ^ 


Exp. 
Time (s) 


050318 


PC 


2.8 


3280 


23431 


050525 


PC 


9.1 


7056 


3560 


050802 


PC 


1.8 


480 


4835 


050824 


PC 


3.6 


6096 


17157 


051111 


PC 


5.0 


5552 


1583 


060418 


WT 


9.2 


448 


1841 


060512 


PC 


1.4 


3680 


58966 



Measured from the BAT trigger time. 

compatible with xspec (version 12.2.1), and response ma- 
trices were taken from version 102 of the UVOT calibration 
files. To create an SED at an instantaneous epoch the count 
rates in these files were set to correspond to the count rate 
of the GRB at the appropriate epoch. These count rates and 
the associated errors were determined from power law model 
fits to the light curve in each filter where the GRB afterglow 
is assumed to decay at the same rate in the UV and optical 
bands. This is justified by the small variations present in the 
decay rate between filters for the GRBs in our sample, which 
have a typical variance of ~ 0.2 and are consistent within 
1(7 errors. Furthermore, the colour evolution expected when 
the cooling break migrates through the optical bands is typ- 
ically observed at later time s than the epochs w e deal with, 
on the order of ~ 10* s (e.g. iBlustin et al.ll200d ). 

Where possible the UV/optical decay index was taken 
from the literature where UVOT data were used in the anal- 
ysis. Otherwise, this analysis was carried out by ourselves, 
with the exception of GRB 050824 and GRB 051111, for 
which there was not sufficient UVOT data to constrain the 
fits to the light curves and, therefore, reported indices from 
J?-band light curve fits were used. For the analysis done on 
UVOT data, both our own and that reported in the liter- 
ature, the decay index was determined from the combined 
UVOT light curve, which was produced by normalising each 




''X (/xm ') 



Figure 1. Galactic (solid), LMC (dashed) and S MC (d otted) 
extinction curves. Parameterisations taken from iPeil lll992l ) 



filter to the 1^-band light curve. The count rate at an instan- 
taneous epoch for each filter was then determined from the 
best-fit temporal decay model using the appropriate normal- 
isation. The light curves were modelled as a power law, with 
the exception of GRB 050525, which required further com- 
ponents to describe the more complex early time temporal 
behaviour l|Blustin et alJ l 20061. The decay rates determined 
for each burst are listed in Table [TJ as well as the references 
used where appropriate. 



2.2 X-Ray Data 

The XRT covers the 0.2 - 10 keV energy range, and all data 
were reduced using the xrtpipeline tool (version 0.9.9). In 
most cases the data used were taken in photon counting 
(PC) mode with the exception of GRB 060418, for which 
the majority of the data take n at the time of interest are in 
window timing (WT) mode (iHiU et all 120041 '). For PC data 
source counts were extracted from a circular region centred 
on the source with an outer radius ranging from 50" to 95". 
In the case of GRB 050802 the data suffered from pile-up 
and, therefore, an annular extraction region was used to ex- 
clude those pixels that were piled-up. The inner radius to 
this was 9.44" (4 pixels) and the outer radius was 106.2" 
(45 pixels). The background count rate was estimated from 
a circular, source-free area in the field of view (FOV) with a 
radius of 118" (50 pixels). For WT mode data, the extrac- 
tion regions used for the source and background count rates 
were 94" slits positioned over the source and in a source free 
region of the FOV, respectively, xselect (version 2.4) was 
used to extract spectral files from the event data in the en- 
ergy ranges 0.3 - 10 keV for PC and WT mode data, which 
is the recommended band to use for compatibility with the 
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current calibration files ^ . Corresponding effective area files 
were created using the xrtmkarf tool (version 0.5.1) and 
binned to at least 20 counts per energy bin. Response ma- 
trices from version 8 of the XRT calibration files were used 
for both WT and PC mode data. 

A large proportion of GRBs observed by XRT have 
shown spectral evolution in the first few thousand seconds of 
emission, during X-ray flares and after early time temporal 
breaks, which could be the transition from internal shock 
to external shock domin ated emission (e.g iBurrows et al.l 
l2005bl : IZhang et all 120061 ). AU X-ray data were taken from 
time intervals where spectral evolution was no longer ob- 
served, and the spectra were normalised to the epoch cor- 
responding to the SED by using the best-fit model to the 
X-ray light curve, in the same way as for the UVOT data. 
The data mode and time intervals used are listed in Table [S] 



3 THE MODEL 

To model the afterglow spectral continuum we tried both 
a power law and broken power law fit, where in the latter 
the change in spectral slope was fixed to A/3 = 0.5 to corre- 
spond to the change in slope caused by a cooling break. In 
addition to this a constraint was also imposed on the break 
energy such that it was within the observing window (i.e. 
0.002 keV < Eb < 10.0 keV). Any fit with a spectral break 
outside this energy range would be equivalent to a power 
law fit to the data. In both the power law and broken power 
law models we included two dust and gas components to 
correspond to the Galactic and the host galaxy photoelec- 
tric absorption and dust extinction. The Galactic column 
density and redde ning in the line of si g ht we re fixed to the 
values taken from IPickev fc LockmanI (|l990h and Schlegel, 
Finkbeiner & Davis (1998), respectively. The second pho- 
toelectric absorption system was set to the redshift of the 
GRB, and the equivalent neutral hydrogen column density 
in the host galaxy was determined assuming solar abun- 
dances. The dependence of dust extinction on wavelength in 
the GRB host galaxy was modelled on the empirical extinc- 
tion laws corresponding to the MW, LMC and SMC. 

The wavelength dependence on dust extinction ob- 
served in these three environments is well reproduced by a 
dust model composed of silicate and graphite grains, where 
variations in the relative abundance and grain size distri- 
bution produce the different extinction laws (jPeil Il992h . 
The parameterisation of the dust extinction laws in the 
MW, LMC and SMC are shown in Fig. [H where Rv = 
Av/E{B -V) = 3.08, 2.93 and 3.16 for th e Gala ctic, SMC 
and LMC extinction laws, respectively 

(lPei|ll99^)- 

The prominence of the 2175 A absorption feature and 
amount of FUV extinction vary between the three curves. 
The MW has the strongest extinction at 2175 A and small- 
est amount of FUV extinction, whereas the SMC has the 
greatest amount of FUV extinction, rising faster than 1/A, 
and an insignificant 2175 A feature. A further difference in 
these three extinction law s is in the amount of reddening 
per H atom (|Drainell2000l ). which is observed to be greatest 
in the MW and least in the SMC. 

^ http: / /heasarc. gsfc.nasa.gov/docs/heasarc/caldb/swift/docs/xrt/ 




Figure 2. Ratio of rest-frame Njj.x to Ay for the three models 
used in the spectral analysis; SMC model (top panel), LMC model 
(middle panel) and Milky Way model (bottom panel). Filled cir- 
cles correspond to power law spectral models and open circles to 
broken power law fits. Dashed and dotted lines show the mean 
Njj X / Av value determined from each spectral model for a power 
law and broken power law fit, respectively. 

We chose these three curves to model the rest frame 
visual extinction in the GRB afterglows to identify the pre- 
dominant extinction properties of the dust within the GRB 
local environment. We refer to each of the spectral models 
as the MW, SMC and LMC model, where the name cor- 
responds to the extinction law used to describe the dust 
extinction properties in the GRB host galaxy. 



4 RESULTS 

The value that we determine in our spectral analysis 
for Nh is an equivalent neutral hydrogen column den- 
sity that results from the amount of soft X-ray absorp- 
tion in the spectrum, where solar abundances are assumed. 
This is dominated primarily b y oxygen K-shell absorption 
([Morrison fc McCammonlll983l ). To distinguish between the 
equivalent neutral hydrogen column density determined 
from the X-ray absorption and the true neutral hydrogen 
column density, we use the notation Nh.x to refer to the 
former. 

The results from our spectral analysis are provided in 
Tableland Table |4] for a power law and a broken power law 
fit, respectively. For five GRBs a broken power law provided 
little improvement to the fit with an F-test probability P 
> 0.08. For GRB 050802 and GRB 060418 a spectral break 
at 1-2 keV and ~ 3 keV, respectively, improved the good- 
ness of fit of the model with an F-test probability ranging 
from 0.003 to 1.7 x 10"^ for GRB 050802 and from 0.02 to 
9 X lO""* for GRB 060418, depending on the model. How- 
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Table 3. Results from simultaneous UV/optical and X-ray power law spectral fits, where only 
UVOT data with rest-frame wavelength A > 1215 Awere used. 



GRB 


Model 


Nh.x " 
10^1 cm"2 


/3 


Av " 


(dof) 


Null Hypothesis 
Probability 


UOUoiO 




'-■^"-0.4,2 


u.yo m u.uo 


U-Oo m U-UD 


iUi yiOj 


u.uo 




LMC 
MW 


1 00+0.53 
^•^o_0 51 
1 go+0-60 


1.04 ±0.04 
n qq+0.06 


0.91 ±0.14 


105 (76) 
154 (76) 


0.02 
0.00 


050525 


SMC 




0.76 ±0.01 


0.16 ±0.03 


29 (33) 


0.67 




LMC 




0.79 ±0.02 


0.23 ±0.04 


33 (33) 


0.48 




MW 




0.78 ±0.03 


0.23 ±0.06 


58 (33) 


0.01 






n qi+O-41 
'-'■''^-0.39 


U.DD m u.uz 


n OA -\- n fi'i 
m U-Uo 


1 0R 






LMC 


1 so+o"" 

1.0U_j) ^2 


0.70 ±0.02 


0.38 ±0.04 


98 (81) 


0.10 




MW 


-'-•="'-0.48 


0.76 ±0.03 


'^■'^0-0.07 


88 (81) 


0.29 


UoUoz4 




1 na+0.42 
J^-U^-0.37 


n n _l_ n no 
U.yo ± U.Uz 


o 1 o _l_ o riA 
0.12 ± 0.04 


O'Z /I \ 

2/ (24j 


0.33 




LMC 


1 14+0.43 
— 0.38 


0.96 ±0.03 


0.18 ±0.05 


27 (24) 


0.29 




MW 


^•^"-0.38 


0.98 ± 0.04 


0.27 ±0.10 


31 (24) 


0.16 


051111 


SMC 
LMC 
MW 


7.71+i-8» 
8.50+i-?9 


1.10 ±0.06 

-I 10+0.09 

1.19 ±0.13 


0.64l»;f3 
0.85^^3! 


13 (11) 

14 (11) 
22 (11) 


0.28 
0.23 
0.03 


060418 


SMC 


0.97tg:g 


0.89 ±0.01 


0.17 ±0.02 


79 (75) 


0.36 




LMC 


1 OQ+0.76 
^•°'-'-0.68 


0.97 ±0.02 


0.38 ± 0.05 


68 (75) 


0.70 




MW 


^•'^-0.71 


0.88 ±0.04 


0.15 ±0.10 


156 (75) 


0.13 


060512 


SMC 


< 0.34 


0.99 ±0.02 




33 (20) 


0.04 




LMC 


< 0.33 


0.98 ±0.02 


0.44 ± 0.05 


37 (20) 


0.01 




MW 


< 0.27 


0.96 ± 0.02 


0.37 ±0.04 


44 (20) 


0.002 



" At the redshift of the GRB 



ever, the column density and the dust extinction at the host 
galaxy do not change significantly between the two contin- 
uum fits, remaining consistent to 2a. 

For the purpose of this paper, where we are interested 
in the properties of the GRB local environment, the factors 
of importance are the dust extinction law that best fits the 
GRB afterglows and the relation between the column den- 
sity and visual-extinction in the GRB local environment. We 
find that the former of these is not affected by the models 
that we use to fit the continuum, and there is no significant 
change in the ratio between the host galaxy column den- 
sity and rest-frame visual extinction. This is illustrated in 
Fig. [2] which plots the value of Nh,x / Av for the best-fit pa- 
rameters determined from a power law and a broken power 
law fit as solid and open circles, respectively. Beyond this 
section we, therefore, refer only to the spectral results from 
the power law fits. In Fig. |3]we show the SEDs and best-fit 
power law models for each GRB. In these figures data points 
at rest-frame wavelengths A < 1215 A, which were not used 
in the fits (see section [2]), are shown as open circles. 

The amount of intrinsic absorption and extinction de- 
termined in our sample varies by over a factor of five in Av 
and by more than an order of magnitude in Nh.x- Typically 
the SMC model results in the smallest Av in the host galaxy 
and the MW model in the largest. This general trend is to 
be expected given the difference in FUV extinction observed 
in the three curves. The MW extinction law is the shallowest 
of the three extinction curves (see Fig. [l| and, in particu- 
lar, has the least amount of dust-absorption in the FUV for 



a given Av. A spectral model with a MW extinction law, 
therefore, requires a larger Av to fit the same data. 

Although the X-ray column density and UV/optical ex- 
tinction are independent components in the fit, the fitted 
value of Nh,x depends on the spectral index, which in turn 
depends on the UV/optical extinction. The best-fit X-ray 
column density will, therefore, vary between extinction mod- 
els. Typically the MW model requires a steeper spectral 
index to compensate for the reduced amount of FUV ex- 
tinction, and this consequently results in larger absorption 
in the soft X-ray. Regardless of these differences the hydro- 
gen equivalent column density determined from our spectral 
modelling is generally consistent at the la level. In the case 
of GRB 050802 the column density is consistent at the 2a 
level between dust models. 

To further investigate the model dependence between 
the amount of X-ray absorption and dust extinction in the 
local environment of the GRB we show the confidence con- 
tours of E{B — V) vs. Nh,x for the best fit power law models 
to each GRB in Fig. 3] For the most part the contours are 
fairly circular, indicating that there is no significant corre- 
lation between Nh.x and Av in our spectral modelling. We 
also tried fitting the X-ray data alone to make sure that the 
UV/optical data were not in any way skewing the column 
densities. We find the best-fit Nh.x value determined from 
our SED spectral analysis to be compatible with the best-fit 
Nh.x value from spectral analysis on the X-ray data alone, 
as shown in Fig. [5] This shows the robustness of our method. 
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Table 4. Simultaneous UV/optical and X-ray broken power law spectral fits. Only using UVOT 
data with rest-frame wavelength A > 1215 A. 



GRB 


Model 


Nh.x " 
10^1 cm~2 


Pi 


(kcV) 


/32 




(dof) 


Null Hypothesis 
Probability 


050318 


SMC 
LMC 
MW 


1 70+0-'*2 

r, r,o+0.41 

^■^o_0 44 
-'-■»y_0.48 


47+0-03 
'-'-^'-0.02 

0.54 ± 0.04 
0.50 ±0.05 


0.003 ± 0.0003 
0.002l0;00i 
< 0.0004 


97+0.04 

1.04 ±0.04 
1.00 ±0.05 


o.56l°:°^ 

0.83lH^ 

91+0-13 


99 (75) 
105 (75) 
154 (75) 


0.03 
0.01 

2e-^ 


050525 


SMC 
LMC 
MW 


1 Qq+0-48 

92+0-30 
^■^^_0 37 

-^■■"-'-0.40 


42+0-" 
'-'-^^-0.10 

30+0-01 

67+0-03 
^-°'-0.04 


029+0-104 

U.U^»_Q 022 

004+°""^ 

U.UU'i_Q 001 

1.048 ± 0.290 


0.9210- 

0.80l°gi 
1 17+0.03 

l-l'-0.04 


0.26 ± 0.04 
0.38lO;«? 
0.24 ±0.07 


27 (32) 
32 (32) 
80 (32) 


0.73 
0.47 

5e-6 


050802 


SMC 
LMC 
MW 


1 14+0-44 
— 43 
-, 70-1-0.51 
^- '•-'-0.49 


0.61 ±0.02 

0.64 ±0.03 
n 79-1-0.03 

^•'^-0.05 


2 970+0-286 
z.»lu_o 34g 

2 990+0-51'' 
z.yyu_o 24g 

3 9fi0+0-538 
o.aou_o 900 


1.11 ±0.02 
1.14 ±0.03 

1 99+0.03 
1-^^-0.05 


0.18 ± 0.03 

0.28in^ 
0.551°;?? 


84 (80) 
81 (80) 
79 (80) 


0.36 
0.43 
0.53 


050824 


SMC 
LMC 
MW 


1 16+0'"' 
^■^"-0.37 
-, 97-1-0.45 
J^-^'_0.44 


47+0-05 

'-'■^'-0.03 

n 50+0-03 
n 5S+0-03 

'-'-^'^-0.04 


0.003 ± 0.001 
0.003 ± 0.001 
0.003 ± 0.001 


97+0 05 

-0.03 

1 00+"°^ 

±.uu_o 04 

1 03+0 03 


16+0-06 

0-251°:^^ 
0.3810-03 


26 (23) 

27 (23) 
30 (23) 


0.32 
0.28 
0.14 


051111 


SMC 
LMC 
MW 


10.041-/70 

8.50ti:g 
11.451^:^3 


o-82«:ro 

0.68«:« 
95+0-18 


040+0-820 
U.U+U_o 034 

< 0.772 

< 0.883 


-1 09+0.24 
l-'-'^-().20 

i.i8oi»:- 
i-45l^;fo 


0.42 ± 0.13 
64+0-26 
t;o+0--*l 


12 (10) 
14 (10) 
23 (10) 


0.25 
0.17 
0.01 


060418 


SMC 
LMC 
MW 


2 05+0-''* 

1 64+0-93 
^•°^-0.89 


0.85 ±0.01 

0.95 ±0.02 
79+0.02 
'J-'^-o.oi 


1 97Q+0.254 
l-^"'-a.l87 

2 100+0-68* 
Z.1UU_0 513 

1 QA7+0-526 
l-o^if _o 150 


1.35 ±0.01 
1.45 ±0.02 
1 29+0-02 


0.17 ±0.02 
0.37 ±0.05 
< 0.03 


68 (74) 
64 (74) 
145 (74) 


0.68 
0.79 

2e-6 


060512 


SMC 
LMC 
MW 


< 0.47 

< 0.42 

< 0.32 


54+0-05 
0.53 ±0.01 
0.51 ±0.01 


0.003 ± 0.001 

0.003 ± 0.001 
< 0.001 


1 04+0-05 

1.03 ±0.01 
1.010 ±0.01 


n t;o+0.06 

57+0.06 

49+0-04 
^-^^-0.003 


32 (19) 
35 (19) 
44 (19) 


0.04 
0.01 
0.001 



" At the redshift of the GRB 



4.1 GRB Host Extinction Laws 

Six of the seven GRBs in our sample were best fit by the 
SMC or LMC model, with the MW model rejected with 
at least 97% confidence for four of these (GRB 050318, 
050525, 051111 and GRB 060512). For GRB 050824 and 
GRB 060418 the MW model is rejected with 84% and 87%, 
respectively. The relatively small amount of extinguishing 
dust in the circumburst environment of those two latter 
GRBs, indicated by the best fit Ay, is likely to be the cause 
for the smaller distinction between the spectral models. Al- 
though the distinction in the goodness of fit between the 
SMC and LMC models is small, the SMC model provides 
the best fit to five of the GRBs in the sample, and only 
GRB 060418 is best fit by the LMC model; = 68 for 75 
degrees of freedom (dof) compared to = 79 for 75 dof for 
the SMC model. 

The afterglow of GRB 50525 was de t ected in aU 
six lenticular UVOT filters (iBlustin et al.1 l2006h . and 
GRB 050 824 ('Schadv fc Campanal [2005h and GRB 060418 
ISchadv fc Falcone had an afterglow detection in all 

but the bluest UV fiher {UVW2). GRB 051111 did not have 
an af terglow detected in the two blue st filters |Poole et ahl 
I2OO5I ). and in the case of GRB 050318 JStill et al.ll2005l ) and 
GRB 060512 ijPe Pasauale et allliooi ) no afterglow was de- 
tected in any of the UV lenticular filters above the 3(t level. 

At the redsh ift s of GRB 050318 [z = 1.44; 
iBerger fc MulchaevI (|2005l )]. GRB 050824 [z = 0.83; 



( Fvnbo et alll2005bl ')]. GRB 051111 \ z = 1.549: iProchaskal 
(2005')] and GRB 060418 [z = 1.49; (|Dupree et al.1 12006*)]. 
the lack of an optical afterglow detection in the bluest filters 
could either be the result of Lya blanketing or high levels 
of dust extinction blueward of ~ 3800 A. GRB 050318 and 
GRB 051111 have Av values that lie at the higher end of the 
distribution observed in the sample, with best fit parameters 
Av ~ 0.53±0.06 mag and Av ~ 0.39lo }o nif-gi respectively, 
if the SMC model is used. GRB 050824 and GRB 060418, on 
the other hand, have at least half this amount in rest-frame 
extinction if the SMC model is used {Ay = 0.12 ±0.04 and 
Av = 0.17 ± 0.02 mag, respectively). The correlation be- 
tween Av and the number of filters in which the afterglow 
is detected could indicate that it is dust present in the local 
environment of these GRBs that contributes to the observed 
dimness of their UV afterglow. 

GRB 060512 was a low redshift burst [z=0.4428; 
iBloom et all (|2006l )]. eliminating neutral hydrogen absorp- 
tion as the cause for the lack of an UV afterglow detec- 
tion. The best fit is provided by the SMC model [x^ = 
33 for 20 dof) yielding a host galaxy extinction of Ay 
= 0.44lg 

MS m^-Si which is comparable to that observed in 
GRB 050318 and GRB 051111. This, therefore, provides fur- 
ther support to the hypothesis whereby dust in the local 
environment of the GRB blocks a large fraction of the UV 
flux emitted. 

In contrast to the other bursts discussed in this section, 
GRB 050525 had very small amounts of local absorption 
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Figure 3. SEDs for seven GRBs at an instantaneous epoch (see text) with best-fit models for each corresponding dust extinction curve 
shown; SMC (dashed), LMC (sohd) and Galactic (dotted). Open circles are data points at wavelength A < 1215 A in the rest frame, and 
therefore not used in the spectral fitting. The positions of 2175 A and 1215 A in the rest frame are indicated. 



and extinction. However, good qualit y data resulting from 
the proximity [z=0.606; B lustin et al.l (2006)] and brightness 
of this burst in the UV and optical bands constrained well 
the spectral fits, and provided a distinction between them. 
The afterglow SED was best fit by the SMC model (x^ = 29 
for 33 dof). 

GRB 050802 is the only burst in our sample where the 
goodness of the spectral fit is improved with the MW model, 
(X^ = 88 for 81 dof). The SED of GRB 050802 flattens 
out at longer wavelengths (Fig. |3]) and this is well fit by a 
model with a dust extinction curve that contains the 2175 A 



absorption feature (i.e. LMC and MW models). The of 
the LMC model fit is still acceptable, with = 98 for 81 
dof, although it is rejected at the 90% confidence level in 
contrast to the MW model, which is only rejected at the 
71% confidence level. 
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Figure 4. Confidence contours for E(B — V) vs. Njj.x from the spectral model fits to the SEDs of the GRBs in the sample. For each 
GRB, the confidence contours shown are taken from the dust extinction model that provided the best fit. The contours are drawn at 
Ax'^ = 2.3,4.61,9.21, corresponding to 68%, 90% and 99% confidence for two interesting parameters. 



5 DISCUSSION 



absence of this feature in the spectra of GRBs indicate it to 
be rare in GRB host galaxies. 



5.1 The 2175 A Absorption Feature 

The 2175 A absorption feature faUs in the UVOT wave- 
length range for the GRB sample used in this paper, and its 
prominence in the GRB SEDs provides information on the 
graphite content of small grains in the surrounding circum- 
burst material. The origin of the 2175 A feature is most likely 
to be carbonaceous material such as small spherical particles 
of graphite (a ^ 30 nm) , which hav e a strong feature at this 
wavelength and of a similar width (jPraine fc Ledl 19841 ). Its 
strength in the MW extinction curve would require ~ 15 % 
of the solar abund ance in carbo n to be present in small par- 
ticles of this size l|Drainell2003l ). whereas its absence in the 
SMC extinction curve can be explained by a differenc e in 
the r elative abundances of graphite and silicate grains (|Peil 
Il992t ). 

The evidence for the 2175 A absorption dip in the 
SED of GRB 050802 suggests a larger abundance of small 
carbonaceous grains in the surrounding environment of 
this burst than is the case for the other bursts. A few 
other GRBs have a lso shown evidence for such a featur e 
[e.g. GRB 970 508: llStratta et all 12004 IXann et all [200^ . 
GRB 991216; l|Vreeswiik et al.l l2006h ] ahhough the usual 



5.2 X-ray Absorption vs. UV/Optical Extinction 

The amount of dust extinction observed in the afterglow of 
GRBs is a measure of the column density of dust grains re- 
sponsible for the absorption of UV and optical photons, and 
the ratio between the X-ray column density and extinction 
gives an estimate of the gas-to-dust ratio in the surround- 
ing environment of the GRB. In Fig. we show the region 
of Av and Nh,x parameter space occupied by our sample 
of bright GRBs for the best fit Av and Nh,x values de- 
termined from the spectral analysis when fitting the SMC 
(top panel), the LMC (middle panel), and the MW models 
(bottom panel) to the data. 

The dashed lines in the figure are plotted as a point of 
reference and correspond to the empirical relation observed 
between Ay and Nh,x in the SMC (top), LMC (middle) and 
MW (bottom). These are determined from the Nh / Ay val- 
ues reported in the literature for each of these environments, 
which are then converted to an Nh,x/Av ratio relating to 
the column density that would be measured from X-ray ob- 
servations of the galaxy if solar abundances were assumed. 
The parameterisation of these lines differ in each panel, and 
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Figure 5. Best-fit Njj,x from SED spectral analysis vs. best-fit 
X from spectral analysis to X-ray data alone. Top, middle 
and bottom panel correspond to the results from the SED spec- 
tral analysis using an SMC, LMC and MW model, respectively. 
In each panel the dashed line corresponds to NH,xO^-'':^y) = 
^-H,x(SED). All data points lie on or very close to this line, il- 
lustrating the robustness of our spectral analysis. 
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8 
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p-(MW) = 0.18 X 10^2 cm" 
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where the factors of | and | accoun t for the l ower metalhci- 
ties observed in the SMC and LMC dPeill 19921). The Nh /Av 
relations are taken f rom I Weingartner fc Draind l|200d ) and 
iPredehl fc SchmittI (ll995D for the SMC and M W, respec- 
tively, and the average of th e ratios found by iKoornneei 
l| 19821 ) and lFitzpatricM j 19851 ) are used for the LMC. 

The data points in Fig.[6]primarily lie to the right of the 
lines of constant Nh,x / Av ■ However, they are confined to 
an area of the Nh,x-Av parameter space much closer to the 
dust-to-gas ratio observed in the MW and Magellanic Clouds 
than the region of space occupied by previous data. This 
is illustrated in Fig. [BJs, where we include the host galaxy 
Nh.x and A v for t wo pre- Swift GRB s amples taken from 
IStratta et all (|2004l ) (SFA sample) and iKann et al.l (|2006l ) 
(KKZ sample), which are shown as open triangles and open 
squares, respectively. Despite apparent differences between 
the Swift sample and previous samples, the distribution in 
Av a nd Nh.x between sa mples remain consistent within 
errors. IStratta et al.l (|2004l ) did not use the LMC extinction 
law to model the afterglow spectrum, and consequently there 
are no triangles shown in the centre panel of Fig. |6}3. 

The X-ray data in the SFA sample are from BeppoSAX 



and were taken hours to days after the prompt outburst, 
with the earliest observation beginning at T + 4 hours. In 
contrast to this, the GRBs in our sample typically have X- 
ray data starting hundreds of seconds after the BAT trigger, 
and the longest delay between the GRB prompt emission 
and the first X-ray o bservations is 1.7 hours (GRB 050824; 
ICampana et al.ll2005l ). The larger signal-to-noise provided at 
early times when the afterglow is significantly brighter im- 
proves the accuracy of our spectral analysis a nd reduces the 
errors associated with the best fit parameters. IStratta et al.l 
(l 2004h point out that the quality of the X-ray data for 
the majority of the bursts in their sample does not allow 
for significant detections of host galaxy absorption in addi- 
tion to Galactic, and in only the cases of GRB 990123 and 
GRB 010222 were the presence of excess absorption robustly 
detected. This is indicated in Fig.[6}3, where the Nh.x value 
for all but two of the SFA sample are upper limits. 

iKann et al.l l|2006l ) focused primarily on optical and NIR 
data, where X-ray column densities used in their analysis 
were taken from the literature. In their sample the earliest 
X-ray observation was still only ~ T -I- 4 hours and the av- 
erage delay was nearly 50 hours from the time of initial out- 
burst. However, for most of their sample the column density 
was determined from higher quality X-ray data taken with 
XMM-Newton or Chandra. For this sample the host galaxy 
Nh.x ranged from an undetectable amount (GRB 021004; 
Moiler et al] l2002l) to (12l^) x 10^^ cm"^ (GRB 010222; 



Stratta et alll2004l ). and the mean is (4.6^1;^) x 10^^ cm' 



which is consistent within 2cr of our sample, which is 
(2.2 ± 0.3) x 10^^ cm"2. Their optical and NIR analysis of 
19 GRBs provided a distribution in Av that ranged from a 
negligible amount up to Av= 0.80 ± 0.29 mag, which is in 
good agreement with our results, for which Av ranges from 
Av = 0.12 ±0.04 to Av = 0.65l^:J^? mag. The mean extinc- 
tion is 0.21 ± 0.04 mags and 0.38 ± 0.02 mags for the KKZ 
and our sample, respectively. 

The high quality of the data in our sample and the 
simultaneous spectral fitting of the X-ray and UV/optical 
data is providing greater constraints on the spectral mod- 
elling, and consequently reducing the errors on the best-fit 
parameters. The results from this are that the gas-to-dust 
ratios in the local environment of GRBs are in better agree- 
ment with those observed in the Milky Way and Magellanic 
Clouds than previous data suggest. The relatively large gas- 
to-dust ratios in GRB local environments indicated by pre- 
vious data were interpreted as evidence of dust destruction 
by the GRB, which would cause the value of Av to decrease. 
However, the results from the Swift data analysis show little 
evidence of this. 

Moreover, when interpreting the Nh.x to Av ratio in 
the GRB local environment compared with that of other en- 
vironments, it is necessary to consider more than just the 
effect of the GRB emission on Av- Photoionisation of the 
gas in the surrounding environment by the GRB X-ray ra- 
diation causes Nh.x to decrease with time, and the extent 
to which the GRB reduces the Av and Nh.x will depend on 
the properties of the GRB and its local environment, such 
as the prompt and afterglow spectral and temporal indices, 
and the density and density profile of the absorbing and 
extinguishing material. 

The lack of evidence for any colour evolution in our 
UV/optical data (see section \2.1\ . provides an upper limit 
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Figure 6. Host galaxy Ay vs. Nh,x- In both plots solid circles represent the host galaxy Ay and Nh,x determined from our spectral 
analysis using the SMC model (top panel), the LMC model (middle panel) and the MW model (bottom panel). The dashed curves are 
the Nu x/-^v ratios for each corresponding environment. This was determined using the Njj / Ay ratios reported in the literature, where 
Njj is converted to the X-ray equivalent Njj x value assuming a metallicity 1/8 and 1/3 solar for the SMC and LMC, respectively (see 
text for details). In the right panel the open triangles and open squares are GRBs taken from the SFA and KKZ sample, respectively. 



on the time by which the GRB no longer destroyed signif- 
icant amounts of dust. i.e. the time at which colour infor- 
mation is available with UVOT, which is typically within 
~ 10'^ s. A limit on the time after which the GRB no longer 
photoionises the surrounding environment can also be de- 
termined by investigating the change in Nh,x over time 
during the early stages of the X-ray afterglow. For this pur- 
pose, we applied spectral analysis on the early time X-ray 
data of the five GRBs in our sample for which there were 
data within the first few 100 seconds of the BAT trigger. 
Of these GRBs only GRB 060418 shows evidence for evo- 
lution in the column density at greater than the 3a level. 
The column density measured from T+84 s to T-)-114 s 
was (1.6 ± 0.10) X 10^^ cm"^, whereas beyond T-l-400 s 
this was i3.9ti l) x 10^^ cm'^. GRB 06041 8 had a very 
large X-ray flare that peaked at T-l-135 s l| Falcone et al.l 
I2OO6I ) and increased the X-ray flux by about an order of 
magnitude in ~ 15 s, which could have caused photoioni- 
sation of the circumburst environment out to greater radii 
from the source. However, it is possible to misinterpret in- 
trinsi c spectral evolution as a change in the column den- 
sity iButler &: Kocevskil I2OO6I ) . The spectral index indeed 
changed from F = 2.75 ± 0.05 to F = 2.15 ± 0.08 between 
the two spectral epochs analysed, and the apparently larger 
Nh,x could be due to intrinsic curvature of the flare spec- 
trum. The column density measured in GRB 060418 no 
longer evolves beyond T-l-300 s, by which time the flare is 
over and, therefore, we determine a typical upper limit of 
a few hundred seconds on the time interval over which the 
GRB photoionises its surrounding environment. 

iPerna fc Lazzatil (|2002l 'l simulated the effect of the GRB 
X-ray and UV emission on Nh.x and Ay over time, where 



they assumed the same initial hydrogen column density 
Nh = 10^^ cm~^ and optical extinction Ay = 4.5 mag 
in all cases, but varied the compactness of the absorbing 
medium and, therefore, also the number density, uh- Their 
simulations indicate that the intense radiation emitted by a 
GRB is capable of photoionising and destroying all gas and 
dust out to a radius of ~ 3 pc within a few tens of seconds. 
Further evolution in the gas and dust column resulting from 
this is, therefore, not expected more than a few hundreds of 
seconds after the peak of the emission, consistent with our 
observations. 

This, therefore, places a lower limit of a few parsecs 
on the scale of the absorbing and extinguishing systems 
detected at the redshifts of GRBs for six of the seven in 
our sample. If the dust and gas were any closer it would 
have been fully d estroy ed and photoionised. Furthermore, 
IPerna fc Lazzatil (l2002l )'s work shows that if the dust and 
gas extends out to a few tens of parsecs the Ay and Nh,x 
measured local to the GRB do not change signiflcantly. This 
would suggest that the host galaxy dust and gas systems 
probed by our sample lie a few tens of parsecs from the 
source, which is consistent with findings by Prochaska, Chen 
& Bloom (2006), who placed a lower limit of J> 50 pc on 
the locations of host galaxy absorption systems from GRBs. 
This suggests that the dust-to-gas ratios measured are ef- 
fectively unaltered by the GRB and a fair representation of 
their local environment. It is, therefore, not so surprising 
that there is no significant deviation in the GRB local envi- 
ronment Nh,x/Av ratio when compared to the Milky Way 
or Magellanic Clouds. 

To verify that the effect of the GRB on the dust and 
gas in its local environment is of little significance, we in- 
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Figure 7. Host galaxy visual extinction against isotropic equiva- 
lent energy in 7-rays for a sample of 16 bursts. Solid data points 
correspond to the GRBs studied in this paper. Open data points 
are GRBs from a pre-Swift sample o f GRBs where the Eigo values 
are taken from [ Bloom et al.l 1I2OO1I') and the Ay val ues are from 
iKann et all ||2006| ') fsguare) and lStratta et al.l ||2004) (triangle). 



vestigate further alternative methods of indirectly detecting 
the process of dust destruction. 



5.3 Alternative Indicators Of Dust Destruction 

Dust destruction, either by sublimation or shattering is de- 
pendent on both the energy of the photons absorbed and 
on the radiation flux, where a less-luminous burst will be 
less effective at destroying the smaller, UV-absorbing dust 
grains. If the radiation emitted by a GRB destroys signifl- 
cant levels of dust in its surrounding environment a correla- 
tion should exist between the energy released by the burst 
and the amount of dust destroyed, which consequently af- 
fects the amount of extinction observed in the GRB after- 
glow spectrum. 

To test this correlation we compute the isotropic equiv- 
alent energy emitted in the 15 - 150 keV 7-ray energy band 
(Eiaoij)) for our sample of GRBs. A k-correction is applied 
using the method described by Bloom, Frail & Sari (2001) 
to convert the prompt energy observed to the same comov- 
ing rest frame bandpass (15 - 150 keV). Fig. [7] shows the 
host galaxy Ay against the k-corrected Eisoiy) for the seven 
GRBs in our sample (solid circles), as well as a further fif- 
teen pre- Swift GRBs (open squares or triangles). All the 
GRBs included in Fig. [7] have spectroscopically measured 
redshifts and estimates of the host galaxy Ay, where the 
host g alaxy Ay for the pre- Swift sample ar e taken from 
either iKann et al.l (I2OO6I ) (open squares) or IStratta et al.l 
(|2004 ) (open triangles). A similar k-correction is applied to 
the non-Swift GRB sample to determine -Eiso(7) in the same 
rest-frame energy range as the Swift sample. 



Twelve of the pre-Swift GRBs had no spectral informa- 
tion and, therefor e, the spectral shape had to be assumed. 
iBand et all (Il993l ') found the GRB prompt emission to be 
well modelled by a smoothed, broken, power law with a soft 
and hard spectral index, a and /3, respectively, and spectral 
break Eq. Those GRBs with no spectral information are as- 
sumed to have a Band spectrum, and their Eisoiy) is taken 
from the average of 54 values, each computed with a dif- 
ferent set of {a,l3,Eo). The set of 54 {a,f3,Eo) c ombinations 
are ta ken from fits to 54 BATSE GRBs done bv lBand et al.1 
(1993). 

Based on the Spearman's rank test the correlation in 
Fig. [7] is not significant, with the log of Eiao and Ay show- 
ing a correlation coefficient of only Vs = —0.21 with a null- 
hypothesis pr obability of P=0.60. This is consistent with 
the findings of lNardini et al.1 (120061 ) ■ whose analysis showed 
a sample of 23 pre-Swift GRBs to show no correlation be- 
tween the isotropic 7-ray emitted energies and the opti- 
cal luminosities. This, therefore, suggests that any dust- 
destruction caused by the GRB intense radiation is not sig- 
nificant enough to affect greatly their observed extinction 
properties, as already indicated in Fig. [6] and discussed in 
section 15.21 



5.4 Extension to Dark Bursts 

By imposing the criteria that the selected GRBs have a 
UV/optical afterglow detected by UVOT, a bias is intro- 
duced that favours afterglows with low source-frame extinc- 
tion. The determined range in host-galaxy extinction may, 
therefore, represent the low end of the extinction distribu- 
tion. Given the association of GRBs with star formation, 
it is likely that the lack of an optical counterpart in some 
bursts is due to large amounts of dust in the host galaxy, 
which blocks out the UV/optical afterglow. In this case the 
distribution of host galaxy Ay may extend much further 
than that which is observed in our sample. GRBs with no 
associated optical counterpart are typically referred to as 
'dark', although the various possible causes for t he lack of 
an optical counterpart make th e term ambiguous (|Ro1 et al.l 
I2OO5I : iRoming fc Masonll2006al ). In this paper we shall use 
the term dark to refer to those GRBs that have no after- 
glow detected at the 3o" level above background in any of 
the UVOT filters within one hour of the prompt emission. 

There are seven Swift GRBs that satisfy both our def- 
inition of dark and which have spectroscopically measured 
redshifts, where z < 5. At higher redshifts than this a GRB 
would always appear dark to UVOT due to the redshifting of 
the Lyman break below the UVOT energy band. To estimate 
the amount of dust in the local environment of this sample 
of GRBs, we initially assume the relation between Nh,x 
and Ay to be linear, and thus use the X-ray column den- 
sity, Nh,x, of the GRB host galaxy to trace the amount of 
dust-extinction affecting the GRB UV/optical afterglow. We 
compare this with the dust extinction estimated for a sam- 
ple of 34 optically bright Swift detected GRBs. This sample 
is made up of those GRBs with an UV/optical afterglow de- 
tected by UVOT and a corresponding spectroscopic redshift, 
up until GRB 060912; this includes our original sample of 
seven UVOT bright GRBs listed in table [31 AU data were 
taken from epochs where no spectral evolution is observed, 
and it was reduced in the same way as described in sec- 
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tion 12.21 The intrinsic column density was then determined 
from a single power law fit with two absorbers; one at z = 
with the column density fixed at the Galactic value, and a 
second one at the redshift of the burst with the column den- 
sity and X-ray spectral slope, Px, left as free parameters. 
The resulting distribution in Nh,x within the two popula- 
tions of bursts is shown in Fig. [8l where the solid histogram 
corresponds to the dark bursts, and the UVOT bright bursts 
are represented by the dashed histogram. Fourteen bright 
GRBs had negligible absorption at the host galaxy, repre- 
sented in the smallest bin of Nh,x- The mean logarithm of 
the X-ray column density, over and above the Galactic col- 
umn is < 20.54 cm"^ and 22.2 ± 0.1 cm"^ for the bright 
and dark population of GRBs, respectively, suggesting that 
dark GRBs reside in denser environments. 

Assuming GRB host galaxies to have an Nh,x and Ay 
relation similar to that observed in the SMC (Eqn.[l}, we use 
the measured Nh,x to determine Ay - Doing this we estimate 
a mean visual extinction of {Av(dark)) = 10.1 mags for the 
dark population of bursts shown in Fig. [S] The most ab- 
sorbed GRB in the sample of dark GRBs was GRB 060510B, 
which had a rest-frame column density Nh,x= (3.54j;Q23) x 
lO'^^ cm~'^ and, therefore, an estimated rest-frame visual ex- 
tinction of Av^ 17.7 mags. However, GRB 060510B had 
an R-band detection of R ~ 21 at around 12 minutes af- 
ter the prompt emission, making the inferred extinction 
rather unrealistic, which suggests that the Nh,x/Av ratio 
of 2 X 10^^ cm"'^ is likely too small for the most absorbed 
of GRBs. 

Instead, we use the host galaxy gas-to-dust ratios de- 
termined in our spectral analysis to acquire a mean value 
of N H ,x / Av , where we use the values determined from the 
SMC model. This corresponds to: 



, Nh.x . 
^ Av ' 



(4) 



Using this relation we estimate a range in rest-frame Av of 
0.4 < Av < 5.3, with a mean of (Av) = 3.0 mags. The 
overlap in rest-frame visual extinction in the UVOT dark 
and bright population of GRBs suggests that several factors 
may contribute to the lack of an optical afterglow. However, 
the mean value of (Av) = 3.0 mags is almost eight times 
the mean extinction observed in the bright GRBs listed in 
Table [3] indicating that the local environments of UVOT 
dark GRBs are much dustier than those with UVOT de- 
tected optical counterparts. This could explain, or at least 
contribute, to the lack of an afterglow blueward of 5500 A. 



6 CONCLUSIONS 

In this paper we use the SEDs of seven optically bright 
GRBs, covering the optical, UV and X-ray energy bands, to 
determine the amount of dust extinction and soft X-ray ab- 
sorption present in the local environments of the GRB. We 
use the Milky Way (MW) and the Large and SmaU Mag- 
ellanic Cloud (LMC and SMC) extinction laws to model 
the host galaxy dust extinction dependence on wavelength. 
These show a decreasing prominence in the strength of the 
2175 A feature, and increasing levels of far-ultraviolet ex- 
tinction, respectively. 

All GRBs but one were best-fit by the SMC or LMC ex- 
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Figure 8. Distribution of the X-ray column density at the red- 
shift of the GRB for a sample of GRBs with (dashed), and with- 
out (solid) a UVOT detected afterglow, where all GRBs have a 
spectroscopically measured redshift. The smallest Nh,x bin rep- 
resents those GRBs with no measurable X-ray column density 
above the Galactic absorption. 



tinction law, with the SMC extinction law typically giving a 
better fit. Only in the case of GRB 050802 did a model with 
a Galactic extinction law at the host galaxy provide a better 
fit, although SMC and LMC extinction laws were rejected 
with only 97% and 90% confidence, respectively. This, there- 
fore, suggests that small graphite grains responsible for the 
2175A feature are not predominant in GRB host galaxies. 

For six of the seven GRBs in our sample an absorption 
and extinction system was detected at the redshift of the 
GRB, which must be located at least several parsecs from 
the source to have survived the intense radiation emitted by 
the GRB. However, the gas-to-dust ratios measured in the 
host galaxies of our sample of bursts are lower than previ- 
ously suggested from analysis of pre-Swift GRBs, and consis- 
tent with those observed in the Milky Way and Magellanic 
Clouds. There is, therefore no evidence of dust destruction 
by the GRB in its circumburst environment provided by the 
spectral analysis alone. However, this does not rule out the 
destruction and photoionisation of dust and gas within a few 
parsecs of the GRB, which is expected to take place within 
the first few tens of seconds of the GRB prompt emission. 

The host galaxy visual extinction observed in our sam- 
ple ranges from Ay = 0.12 ± 0.04 to 0.651!^;^^ mags. How- 
ever, the requirement that the GRBs analysed have a de- 
tected optical counterpart introduces selection effects that 
favour GRBs with a low Av- Using the amount of host 
galaxy Nh,x as an indicator of the amount of visual extinc- 
tion undergone in the GRB local environment, we estimate 
that those GRBs with no optical counterpart have, on aver- 
age, values of Av almost a factor of eight larger than those 
observed in the sample of bright GRBs studied in this pa- 
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per. This would suggest that the host galaxies of dark GRBs 
are intrinsically dustier than those of GRBs with UVOT de- 
tected optical counterparts, which could account for the lack 
of an optical afterglow. 
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